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In a recent review by Anderson and coworkers-, it was pointed out that an early resonating valence 
bond (RVB) theory is able to explain a number of unusual properties of high temperature super- 
conducting (SC) Cu-oxides. Here we extend previous calculations-^^- to study more systematically 
the low energy physical properties of the plain vanilla d-wave RVB state, and to compare the results 
with the available experiments. We use a renormalized mean field theory combined with variational 
Monte Carlo and power Lanczos methods to study the RVB state of an extended t — J model in 
a square lattice with parameters suitable for the hole doped Cu-oxides. The physical observable 
quantities we study include the specific heat, the linear residual thermal conductivity, the in-plane 
magnetic penetration depth, the quasiparticle energy at the antinode (ir,0), the superconducting 
energy gap, the quasiparticle spectra and the Drude weights. The traits of nodes (including Jcf, the 
Fermi velocity vf and the velocity along Fermi surface v?), as well as the SC order parameter are 
also studied. Comparisons of the theory and the experiments in cuprates show an overall qualitative 
agreement, especially on their doping dependences. 

PACS numbers: 74.20.Fg, 74.72.-h, 71.10.-w 



1,5 



X 



INTRODUCTION 



Since the discovery of high transition temperature su- 
perconductivity (HTSC) in cuprates in 1986, there have 
been enormous efforts in both experiments and theories 
to understand the mechanism of the superconductivity 
and their unusual physical properties. One of the ear- 
liest theory is the resonating valence bond (RVB) the- 
ory proposed by Anderson 2 . In that theory, the op- 
erative element in the electronic structure of this class 
of compounds is the square planar CuC>2 lattice. The 
parent compound such as La 2 Cu04, where the Cu is 
stoichiometrically bivalent Cu 2+ with one hole per Cu 
site, is an antiferromagnetically coupled Mott insula- 
tor. Chemical doping such as the partial replacement 
of La by Sr introduces additional holes on the CuOi 
plane. The minimum microscopic model to describe 
the low energy physics has been proposed to be the 
t — J model or its variant extended t — J model, which 
includes an antiferromagnetic spin coupling and a ki- 
netic energy term for the hole motion-^. Anderson pro- 
posed a doped spin liquid of spin singlets, or the bond 
spin singlet resonating between many configurations^ 
This concept explains many unusual properties of the 
cuprates, as emphasized in a recent review by Ander- 
son and coworkers^. More quantitatively, in the sim- 
plest RVB theory, namely its plain vanilla version, the 
RVB state in the cuprate is described by a Gutzwillcr 
projected d-wave BCS wavefunction, whose parameters 
are determined variationally either by using a renormal- 



ized mean field theory (RMFT)£ or by variational Monte 
Carlo method (VMC) numericallji£i&, or by other field 
theory methods-. Recently the Gutzwiller RVB wave- 
function approach was applied to the strong coupling 
Hubbard model by Paramekanti, Randeria and Trivedi^, 
who used careful numerical methods to calculate several 
quantities of direct experimental relevance. Both results 
for the Hubbard and t — J models turn out to corre- 
spond well to some experimental phenomena observed in 
cuprates. The plain vanilla RVB theory has recently been 
extended to study the scanning tunneling microscopy, the 
angle resolved photoemission spectroscopy (ARPES) 10 
and the Gossamer superconductivity— in the Hubbard- 
like models at the half electron filling. In view of the 
preliminary success of the plain vanilla RVB theory, it 
is desirable to extend previous calculations and to ana- 
lyze more experimental data so that a more systematic 
and comprehensive comparison between the theories^iii 
and the experiments^* can be made on more observable 
quantities. 

In the present paper, we extend the previous works 
of Zhang et al£ and of Paramekanti et ali^ to carry out 
more systematic calculations on the low energy physical 
properties of the plain vanilla d-wave RVB state. We use 
a renormalized mean field theory combined with varia- 
tional Monte Carlo and power Lanczos methods (PL)i5ii£ 
to study the RVB state of an extended t — J model in a 
square lattice for parameters suitable for the hole doped 
lanthanum and yttrium Cu-oxides. Our main focus is 
on the microscopic calculations of the key parameters for 
the nodal quasiparticles in the d-wave RVB state, namely 
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the Fermi velocity vf and the velocity along the Fermi 
surface v^- From these quantities, we calculate a number 
of low energy physical properties including the specific 
heat, the linear residual thermal conductivity, and the 
in-plane magnetic penetration depth. We make exten- 
sive comparison between these calculations with a very 
broad spectrum of types of reported experiments, and 
find qualitative agreement, especially on the doping de- 
pendences of these properties. The discrepancy between 
the theory and the experiments is mostly on the abso- 
lute values of these quantities, which may be attributed 
to a factor of 2 ~ 4 times larger in the value of i»2 in the 
theory. We also calculate the quasiparticle energy gap at 
the antinode (n, 0), the SC energy gap, the quasiparticle 
spectra and the Drude weight, and find good agreements 
with the experiments. 

The paper is organized as follows. In section II, we 
describe the microscopic model and the methods we use 
in our calculations. In section III, we calculate the basic 
parameters of the nodal quasiparticles. In section IV, we 
discuss the nodal physics and make comparison of the 
theory and experiments on a number of low energy phys- 
ical properties. In section V, we calculate other physical 
quantities and compare with the experiments. A brief 
summary is given in section VI. 

II. MODEL AND METHODS 
A. model 

We consider an extended t — J model including a near- 
est neighbor (n.n.), a second n.n. and a third n.n. hop- 
ping terms in a square lattice, 

H = H t + Hj (1) 

Hj = J^5i • Sj 

(hi) 

In the above Hamiltonian, a constraint of no double occu- 
pation of electrons on each site is implied: c\ a Ci ytT < 
1. The summations (i, j) , and run over the 

n.n., second n.n. and third n.n. pairs, respectively, i, 
t' and t" are their corresponding hopping integrals re- 
spectively. We choose t and t" to be positive and t' to 
be negative, appropriate for the hole-doped cuprates; J 
is the superexchange coupling between the n.n. spins. 
Table I lists the parameters {t,t' /t,t" /t, J/t} used 
in our calculations for mono-layered La 2 - x Sr x CuO^ 
(LSCO) and bi-layered cuprate YBa 2 Cu 3 7 - x (YBCO) 
or Bi 2 Sr 2 CaCu 2 Os+ x (Bi-2212). These parameters ap- 
pear consistent with the band-structure calculations 17 
and also with the experimental analyses such as the 
topology of large Fermi surface reported in ARPES 1819 , 



the inelastic light scattering 2 -, neutron scatteringSi* 2 ^ 
and two-magnon Raman scattering 24-25 - 2 - experiments. 



TABLE I: Parameters {t,f /t,t" /t, J/t} 
used for La 2 - x Sr x CuOi and Y Ba 2 Cu 3 7 - x 
(BiiSrvCaC^Os+x) in the renormalized mean field 
theory. 





t(eV) 


t'/t 


t"/t 


J/t 


La2- x Sr x CuOi 


0.3 


-0.1 


0.05 


0.3 


YBa 2 Cu 3 7 ^ x 
Bi 2 Sr 2 CaCu 2 O s + x 


0.3 


-0.3 


0.2 


0.3 



We use a variational projected d-wave BCS state or the 
d-RVB state to study the ground state and elementary 
excitations of the models The trial ground state is of the 
form, 

|*l> =Pg\^bcs), (2) 

where the Gutzwiller projection operator Pq = YiiO- ~ 
n^riii) is to ensure the constraint of no double occupa- 
tion of electrons on any lattice site. The BCS state is of 
the standard form, given by 

\^bcs) = Y[(u k + v k c{^_ kl ) |0) 

k 

where |0) is the vacuum, and Uk and Vk are the varia- 
tional parameters satisfying the normalization condition: 

\uk\ + Kl = l- 

In this paper, we use two complementary methods to 
carry out the Gutzwiller projected variational calcula- 
tion. One is the renormalized mean-field theory (RMFT), 
which takes into account of the Gutzwiller projection by 
a set of renormalization factors^ The other is the varia- 
tional Monte Carlo (VMC) method which computes the 
quantities numerically, followed by a further improve- 
ment of the variational wavefunction by using the power 
Lanczos (PL) methodi 5 ^^ to eliminate or to reduce the 
bias in the variational approach. It is well known that 
the variational calculation often overestimates the effect 
of superconductivity in the true ground state, and the 
variational calculation usually leads to a larger A& 

B. Renormalized mean field theory 

The RMFT is a Hartree-Fock like mean field theory 
to approximately treat the projection operator in the 
Hamiltonian JTJ. In the RMFT, we apply the Gutzwiller 
approximation to replace the effect of the projection op- 
erator by a set of renormalization factors, which are de- 
termined by statistical counting 3,27 . The variation of a 
projected state for Hamiltonian H is then approximately 
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mapped onto that of the corresponding unprojected state 
for a renormalized Hamiltoniani^. This method was ini- 
tially developed by Gutzwiller to study possible ferro- 
magnetism in strongly interacting systems^. It was later 
applied by Brinkman and Rice to study the metal insu- 
lator transition 28 , and by Vollhardt to study the Fermi 
liquid theory of Helium-3^. 

Let (Q) be an expectation value of Q in the RVB state 
and (Q)o be an expectation value of Q in the BCS 
state l^scs); then the expectation values of the hopping 
term and the spin-spin correlation in the RVB states can 
be written in terms of those in the BCS state, 



(Si ■ Sj) = g s (Si ■ Sj)o 



(3) 



where gt and g s are the two renormalization factors for 
the kinetic and the spin-spin superexchange terms re- 
spectively, they are given by£ 
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9t 



-, g s = 



1 + 8' as {l + Sy 



with 5 the hole concentration. The evaluation of H in 
the RVB state is then mapped onto the evaluation of the 
renormalized Hamiltonian H' in the corresponding BCS 
state, with H' given by 

H' = g t H t +g s Hj (4) 

The variational energy of the system is then given by 

W = (H) = (H') Q 

In this paper, we shall only consider even parity SC state, 
namely |v-fc| 2 = \vk\ 2 , and v k u k = u*_ k v~ k . We obtain 

W = 2g t Y J W\ 2 e(k) 

k 

+ f? v k-k'(\v k \ 2 \v k >\ 2 + u k vkv* k ,u* k ,) 



k.k' 



where 



e(fc) = — 2t(cos k x + cos k y ) — 4i'(cos k x cos k y ) 

— 2i"(cos 2k x + cos 2k v ), 
3 

Vfe-fe' = -- J[cos(k x - k' x ) + cos(k y - k' y )] 

with V the total number of lattice sites. The total num- 
ber of electron operator N e — J2ka c ka c k<y has an ex- 
pectation value of (N e ) = 2^ fc |vfc| 2 , so that the hole 
concentration 



8 = 1-2^K| 2 /V 



(5) 



Let /i be the chemical potential, the quantity we wish to 
minimize is W — (H — fiN e ) , or 



W' = W- 2^Y^ 



Vk\ 



The variation is carried out respect to v k , u k and 8 for 
fixed /i. Carrying out this variational procedure, we find 
that 



\Uk\ 



(l-£(k)/E(k)), 
{l + £(k)/E(k)), 



u k v k = A(k)/2E(k) 



(6) 



with 



£7(fc) = Vr(fc)+|A(*)r 



The parameters £(fe) and A(fc) are related to the particle- 
hole and particle-particle pairing amplitudes which are 
introduced below in Eci. |7l8j . E(k) turns out to be the 
energy of a quasiparticle in the SC stated We define 



Xr = J2( c L C *+t*) 



(7) 
(8) 



with r = x, y, the n.n. unit vector. For the d x 2_ y 2 
pairing symmetry, A x = -A„ = A , X x = X y = Xo> and 
£(fe) and A(fc) have the forms 

£(fc) = g t e(k) - Jl - x(cosk x + cosk y ), (9) 
A(fe) = A(cos k x — cos k y ) (10) 

where A = (3.g s J/4)A , x = (3g s J/4)x and ju = 
/it + d (H') /NdS. The mean fields Ao and Xo can be 
determined by solving these self-consistent equations jSJ- 
ITU] . The SC order parameter is defined as 



Asc(-Ry) = (4i c ]l ~ c 1| c Jt) 

which is related to the variational parameter Ao in the 
Gutzwiller approximation^^ 



A S c = gt Ao 



Variational Monte Carlo method 



(11) 



In the VMC calculation, we first rewrite the wavefunc- 
tion Eq. in the Hilbert space with fixed number of N e 
electrons doped with even number of n holes, 



*RVB)=PG§2^ct,cl kl )( 
k k 



I 0) (12) 



with 



Wfe 



A M c(fc) 



u k e M c(k) + y/e M c(k) 2 + |A MC (fe)| 2 
ejvfc(fc) = — 2t(cos k x + cos k y ) — At v cos k x cos k y 

— 2i"(cos2fc x + cos2fc a ) — fx v , 
Ajvfc(fc) = 2A v (cosk x — cosky) 
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where A„ and fi v are variational parameters, with A„ 
related to the d-wave SC order parameter and [i v similar 
to the chemical potential. Note that we have used sub- 
script MC to distinguish the parameters here from those 
adopted in the section B, and that we have included two 
additional variational parameters t' v and t", which are 
usually not equal to the bare values t' and t" because the 
constraint strongly renormalizes the hopping amplitude. 
That is to say, the form of eMc(k) in the variational 
wavefunction can be different from the dispersion func- 
tion of the non-interacting electrons. These variational 
parameters determine the Fermi surface topology. Then, 
the quasi-particlc excitations are created by adding holes 
into Eq.(H2J): 

I = Pgc^J: -clA^-")/ 2 - 1 | 0). 

k Uk 

(13) 

From Eq. l|13|) we calculate the energy dispersion for a 
given doping density by using VMC. The system used 
in this paper is of 12 x 12 sites with periodic bound- 
ary conditions. 1 We the n fit the quasi-particle en- 
ergy with the formula a ■ — b to determine 
the renormalized parameters, with a and b the fitting 
parameters. Additionally, in order to eliminate the bias 
introduced in the trial wavefunction method, the power- 
Lanczos methocU£*A£ which is a hybrid of the power and 
the variational Lanczos method is used to further improve 
the trial function. In the power method it can be easily 
shown that if a trial wave function is not orthogo- 
nal to the ground state, (W — H) m |\1>) is proportional 
to the ground state wavefunction as the power m ap- 
proaches infinity. W is an appropriately chosen constant 
to make the ground state energy the largest eigenvalue 
of the W—H matrix. In our calculation, the first order 
Lanczos method, i.e., m = 1 is used and the improved 
trial wavefunction is \PL1) = (1+ C\H) l^). G\ is a new 
variational parameter. The results described below de- 
noted as PL1 are calculated with the trial wavefunction 

(1 + dfOI*). 



III. BASIC PARAMETERS 

In this section, we discuss the parameters of nodal 
quasiparticles in the d-wave SC Cu-oxides and make com- 
parisons between the theory and experiments. It has 
been well established in experiments that the cuprate 
superconductivity has d x i _ y 2-wave pairing 3 ^. There are 
four nodes in the fc-space, where the quasiparticle dis- 
persion Ek = \J (,k + A^ approaches zero. The typical 
Fermi surface (FS) in HTSC is shown in Fig(I] together 



The VMC results reported here are for H in Eq.0 with an 
additional term of (— l/4)jiinj in Hj. This additional term is a 
constant at the half-filled. 




k 

x 

FIG. 1: (Color online) Illustration of the Fermi surface for 
LSCO (dashed square) and for YBCO (solid line), and the 
location of the gap nodes fcp(±l,±l). The "Fermi velocity" 
up and the "gap velcoity" V2 are defined as the slopes of 
the quasiparticle energy along and perpendicular to the nodal 
direction. Vp and V2 specify the Dirac cone for the nodal 
quasiparticle dispersion. 

with the "Fermi velocity" Vp and the "gap velocity" v 2 
which are defined as the slopes of E^ along the direc- 
tion perpendicular and tangential to the Fermi surface 
at the nodes kp(±l, ±1). \J2kp is the Fermi wave-vector 
along the diagonal direction. Different from the conven- 
tional s-wave pairing symmetry, low energy quasiparti- 
cles in the vicinity of these nodes can be easily excited 
by thermal fluctuation, impurity scattering, or disorder 
effects. These low energy nodal quasiparticles predomi- 
nate physical properties of HTSC at low temperatures. 

A. Fermi wave-vector 

From the ARPES data, it is known that a transi- 
tion from a hole-like Fermi surface centered at (tt, tt) to 
an electron-like Fermi surface centered at (0,0) occurs 
slightly above the optimal doping in both LSCO and Bi- 
2212, meanwhile the Fermi wave-vector kp shifts just a 
little biti 18 i 31 i 32 For Bi-2212, ARPES experiments 3 ^ sug- 
gest that Uf is weakly doping-dependent and \/2kp — 
0.434-Ui For optimally doped YBC0 6 . 9 5, V^kp ~ 
O.biA -1 ^^ and for underdoped La 2 -xSr x Cu 2 OA (x = 
0.063), V2k F = O^i- 1 ^ with the lattice constant 
a = 3.8A. These experimental data are shown in Fig|3 
compared with our theoretical calculation where kp is de- 
termined by £(fc) = along the diagonal direction. For 
all the method considered, we found that fcj? decreases 
with increasing doping and kp decreases more rapidly in 
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VMC: O 
PL1: 

ARPES: ★ (ref. 37) 
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(b): YBCO (Bi-2212) 
ARPES: 
■ (Bi-2212 ref. 33,34) 
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0.1 



0.2 



0.1 



0.2 



3.0 
2.5 
< 2-0 

I 

1.5 

> 

1.0 
0.5 
0.0 



(a) : LSCO 




(b): YBCO (Bi-2212) 




★ 




*★ * * 




i ^ ft ^ ^ 












RMFT: 






VMC: O 






PL1: 




Variational theory: O (ref. 4) 


ARPES: ■*• (ref. 38) 




ARPES: ★ (ref. 38) 


0.1 0.2 




0.1 0.2 


5 




5 



FIG. 2: (Color online) Comparison of the doping dependence 
of the Fermi wave-vector kp obtained from our theoretical 
calculation with those obtained by the ARPE S??'? 4 '??'??'?? for 
(a) LSCO and (b) YBCO (Bi-2212). Theoretical results are 
obtained using parameters listed in Table I for LSCO and 
YBCO (Bi-2212) with the in-plane lattice constant a = 3.8A 
RMFT: renormalized mean field theory; VMC: variational 
Monte Carlo; PL1: Power Lanczos to the first order. Note 
that, the values from VMC and PL1 may be very close, and 
get overlapped with each other in some plots following. 

YBCO (Bi-2212) than in LSCO. The values of k F for 
the underdoped LSCO and optimally doped Bi-2212 and 
YBCO agree qualitatively with the experimental data. 



FIG. 3: (Color online) Fermi velocity Vf vs hole concentration 
5. The results from variational theory in Ref^ were obtained 
with the correction of order of O(Jft) included for the Hub- 
bard model. (leV - A ~ 152km/s) 



the same result. In the optimally-doped and overdoped 
regimes, this trend does not deviate greatly from the ex- 
perimentally observed universality. However, the value 
of vf appears underestimated in the RMFT for the ex- 
tended t-J model compared with the experimental data 
and with that obtained for the Hubbard model by includ- 
ing the correction of order of O(Jft) reported previously 
by Paramekanti et alii 



B. Fermi velocity vf 



C. "Gap velocity" vi 



In the vicinity of the gap nodes, the quasiparticle dis- 
persion can be expressed as 



E(k) = 



"II 



where vf (1^2) and ku (fcj_) are the components of the 
velocity and wave-vector perpendicular (parallel) to the 
Fermi surface, respectively. 

The Fermi velocity extracted from the slope of the 
quasiparticle dispersion obtained by the ARPES is found 
remarkably universal, independent of the doping concen- 
tration, vp ~ 250 ~ 270fcm/s within an experimental 
error of 10-20%,2iS. 

In the RMFT, vp and v-i are given by the following 
equations 



= V2\sm(k F ) 



2(5t* , , 
V2 = V2Asm(kp) 



^x)-4g t \t'\cos(k F ) + 8g t t' 



cos 



(14) 
(15) 



Shown in Figgis the value of vf obtained from En.|14|. 
Vp increases with doping. The VMC gives essentially 



The "gap velocity" v-i is the slope of the SC energy 
gap along the Fermi surface at the gap node. Together 
with kp and vp , V2 specifies the Dirac cone for the nodal 
quasiparticle dispersion. Among them u 2 plays a cru- 
cial role in determining the nodal physics of HTSC. This 
is because hp and vf are rather universal, depending 
weakly on the doping concentration. Furthermore, V2 is 
much smaller than vf, and many physical properties are 
related to V2 in the form of the ratio vf/v2, therefore 
a small variation in V2 may lead to a drastic change of 
vf/i>2, hence of some physical quantities. 

Experimentally, V2 is difficult to be determined accu- 
rately. It depends strongly on the doping concentra- 
tion and other material properties. A number of exper- 
iments may be used to extract V2- These experiments 
include ARPES 18 - 3 i, the temperature dependence of in- 
plane magnetic penetration depth A(T)iS, the electronic 
specific heat C e ^2^i^2^, and the linear residual ther- 
mal conductivity «q/T|t=u 36 ' 44145 ' 46 - The linear residual 
thermal conductivity is robust against renormalization 
due to quasiparticle interactions and vertex corrections. 
In the SC state, kq/T\t=o vf/v2 + V2/vf is universal 4 ^ 
and it does not depend on the impurity scattering rate. 

Fig0| shows the doping dependence of V2 and the ratio 
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FIG. 4: (Color online) Panels (a,b) : "Gap velocity" «2 
vs hole concentration 5. For optimal doped YBCO (Bi- 
2212), «2 = 10 ~ 20km/s was reported by various kinds of 
experiments^^Siik. The experimental data indicated in panel 
(a) is achieved by measuring the magnetic field dependence 
of the specific heat on LSCO at the zero temperature limite— . 
Panels (c,d): the ratio vf/v2 vs hole concentration S. The 
doping dependence of the ratio vf/v2 is similar to that ob- 
served in the thermal conductivity experiments (see FigEJ. 



vf/v2- The calculation shows that V2 drops quasi- linearly 
with increasing doping. This behavior is consistent with 
the recent experiment of the magnetic field dependence 
of the specific heat on LSCO 43 . Our calculation gives 
v 2 = 20 ~ 30km/ s at an optimal doping, which is larger 
than the experimentally reported value of 10 ~ 20km/ s 
for YBCO (Bi-2212) obtained from the thermal conduc- 
tivity measuremenlj 4 ^^ and ARPES'Mi 3 ^, and of 7km/ s 
for LSCO obtained from the measurements of magnetic 
field dependence of the specific heat at the zero temper- 
ature limiti 3 .. In the RMFT, v F /v 2 = 1 at zero doping. 
As doping increases, vp /v2 increases rapidly. The dop- 
ing dependence is qualitatively consistent with the data 
reported in the thermal conductivity measurements^^ 

In our theoretical calculation, despite the great differ- 
ence between the hopping integrals t' and t" for YBCO 
(bi-2212) and LSCO, the values of i>2 are only slightly dif- 
ferent. This result is also in qualitative agreement with 
the thermal conductivity measurements. In the next sec- 
tion we will use our theoretical result of V2 and Vf to 
extract some physical observables and compare with ex- 
perimental results. 



IV. NODAL PHYSICS 

In HTSC, d-wave pairing symmetry leads to a dome- 
like quasiparticle dispersion around the gap nodes. In 
the SC state, the gapless quasiparticle excitations in 
the vicinity of nodes dominate low temperature physical 
properties. It is of fundamental importance to explore 
physical properties of these quasiparticle excitations. 

In the clean limit, the density of states (DOS), p(u>), 
of low lying quasiparticles near the nodes is linear, 



p(u>) 



1 



7T VpV 2 



(16) 



The linear coefficient of p(to) is inversely proportional to 
the nodal velocities vf and V2- This linearity in energy 
of p(u>) leads to many unconventional physical behaviors 
such as the quadratic electronic specific hea t 34 ' 3 9i 4 9i 4 9 , 
the linear residual thermal conduct ivitj* 3 ^ 4 * 4 ^ 4 ^ and the 
linear decreasing of superfluid densit}i£2*£i. Experimental 
observations of these behaviors have provided some of the 
early evidences for unconventional d x i_ y i pairing sym- 
metry in HTSC. The nature of the interactions of nodal 
quasiparticles is not so clear in HSTCi^ 4 ^ Some have 
used a renormalization factor to describe the effect of 
quasiparticle interactions on the electronic specific heat 
and on the in-plane magnetic penetration depth^aa& In 
this paper, we shall neglect quasiparticle interactions and 
set the renormalization factor to unity. 



A. Electronic specific heat 

The linear low energy DOS p(uj) leads to a quadratic 
temperature dependence of the low temperature elec- 
tronic specific heat in the HTSC, given by 



jT = aT 2 



21.6 k% 



1 



r, 2 



VFV2 



(17) 



FigEl compares our theoretical results with the exper- 
iments for LSCO and YBCOi 39 i 40 i 41 The experimental 
result for LSCO shows a general tendency to increase 
as doping increases, and the rapid increase of a in the 
overdoped region might be due to the Fermi level cross- 
ing of the flat band at (tt, 0), which yields an addi- 
tional channel to thermally excited quasiparticles. In 
our theoretical results, the doping dependence of a is 
similar for LSCO and YBCO. The values of a are about 
0.01 ~ 0.03mJ/Mol-K 3 which are comparable to the ex- 
perimental value of YBCO) 4 ^^ but much smaller than 
the value of LSCQ22. 



B. Thermal conductivity 

In the presence of a small amount of disorder or im- 
purities, the nodal quasiparticles are delocalized and can 
carry both heat and charge. For dilute non-magnetic 
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LSCO 



YBCO (Bi-2212) 



FIG. 5: (Color online) The quadratic coefficient of the elec- 
tronic specific heat: a — C e i/T 2 vs hole concentration S. 



impurities, there will be a residual normal fluid due to 
these delocalized and conductive quasiparticles. The 
most striking property of this conduction mechanism is 
the universal limit, i.e. the quasiparticle transport is in- 
dependent of the scattering rate as T — > 0. With in- 
creasing the impurity concentration, the mean free path 
is reduced, but the normal fluid density increasesi 47 i 52 
In the SC state with a random distribution of impuri- 
ties of an energy scale Ei m < ksTc, the low temperature 
thermal conductivity is linear and is given by 



~T 



T=0 



k 2 B n vf 



V2 



3h d 



(— + —) 



V2 



VF 



(18) 



where d/n, the stacking distance between two nearest 
neighboring CuOi planes, has the values of 6.6^4, 5.8.4., 
1.12A for LSCO, YBCO, Bi-2212, respectively. This for- 
mula is obtained within the self-consistent T-matrix ap- 
proximation, and it may break down if the impurity scat- 
tering gets too strong^ This universal behavior of the 
thermal conductivity provides a robust and direct mea- 
surement of Vf/v2 in the SC state. 

FigEl shows our theoretical results of kq/T\t=o com- 
pared with the experimental results for LSCO and YBCO 
(Bi-22 1 2). 36 i 45 i 46 i 49 i 54 Experimentally, above a critical 
doping S pc , both LSCO and YBCO (Bi-2212) are thermal 
metals and kq/T\t=o increases steadily as 8 in the under- 
doped regime and very rapidly in the overdoped regime. 
Such observation strongly supports the notion that there 
are well defined nodal quasiparticles in the clean limit. 
The difference of the residual thermal conductivity be- 
tween LSCO and YBCO (Bi-2212) is much smaller com- 
pared with the electronic specific heat shown in Fig|SJ In 
the lightly underdoped regime 6 < 5 pc , the low temper- 
ature behavior of Ko/T remains uncleaB 4 *!^. However, 
it is clear that k /T -> as T -> in LSCO^. 6 . The 
thermal insulating behavior in LSCO is probably due to 
the localization of quasiparticles due to disorder effects. 

In our calculation, the SC state and the delocalized 




FIG. 6: (Color online) Linear residual thermal conductivity 
ko/T vs hole concentration S. 



quasiparticles are assumed to prevail in the heavily un- 
derdoped region. The theoretical results deviate from the 
experimental ones by a factor of 2 ~ 4. We attribute this 
discrepancy to the overestimated gap velocity V2 in the 
theory. 



C. In-plane magnetic penetration depth 

The magnetic penetration depth A(T) is related to the 
superfluid density p s by 



Ps (T) p s (0) Pn (T) c 2 



1 



■in- 



47re 2 A 2 



where m* is the effective mass of the charge carriers, 
and will be assumed to be doping independent, and p n 
is the normal fluid density^ At low temperatures, p n 
is contributed from the thermally excited quasiparticles 
near nodes, and can be given by 



Pn(T) 



nr 



21n2\ nvp ksT 



7T J dv2 h 2 



The linear temperature coefficient of p s (T) / m* is pro- 
portional to Vf/v2- 

At low temperatures, the temperature dependence of 
A(T) is very weak, and A(0) is about several thousands 
angstromsiiSi£ii£Lii£i£2iSI The first and second derivative 
of the penetration depth with respect to temperature can 
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be approximately expressed as 
dX(T) 



e 2 kg n vp 



dT 



A J (0)41n2--|--^ 



d \-2| 



= -8 In 2 



C 2 H 2 d V2 

; 2 kg n vf 



c 2 h 2 d 



V-2 



(19) 
(20) 



Panels (a) and (b) of Fig l7l9l show the zero tem- 
perature in-plane magnetic penetration depth A(0) and 
A _2 (0). Experimentally, as S increases, A(0) in LSCO 
monotonically decreases^Si, while A(0) in YBCO and 
in Bi2Sr2Cai- x Y x Cu20s+s increases with doping in 
the overdoped regionSi^. The experimental results of 
A~ 2 (0) in LSCO and underdoped YBCO (Bi-2212) show 
a linear doping dependence, supporting the idea that the 
zero temperature superfluid density p s (0) is proportional 
to the doping concentration in the underdoped region. In 
our RMFT, in the SC phase, A _2 (0) is nearly linear with 
the hole doping and A(0) diverges at zero doping within 
the approximation that all optical spectral weights are 
condensed to the zero energy in the t-J model. Our the- 
oretical results of A(0)(A~ 2 (0)) agree with the experimen- 
tal data for LSCO and YBCO in the underdoped region. 
In Bi2Sr2Cai^ x Y x Cu20s+s, our theoretical results show 
a discrepancy with the experiments. 

Panels (c,d) of Fig l7l9l show the derivatives of the pen- 
etration depth with respect to temperature, d\(T)/dT 
and d\~ 2 /dT. In the underdoped or slightly overdoped 
region, dX(T)/dT decreases with increasing doping in all 
three compoundsi^liiiii&Li&SiSi In the heavily overdoped 
region d\(T)/dT increases with doping in both LSCO 
and Bi2Sr2Cai- x Y x Cu20%j r g. In LSCO, d\~ 2 /dT in- 
creases steadily with dopingiSLSl In YBCO, an oppo- 
site tendency was observed^ in the underdoped region 
S < 0.10. The anomalous increase in the underdoped re- 
gion was previously shown to qualitatively agree with the 
behavior resulted from the d-density wave statei£2*§2i§i 
In our RMFT, similar doping dependence of d\{T)/dT 
is obtained. However, there is a great discrepancy on the 
absolute values between the experiments and our theoret- 
ical results. We argue that some other mechanisms may 
be responsible for the large value of dX(T)/dT observed 
in experiments 
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V. OTHER PHYSICAL PROPERTIES 



FIG. 7: (Color online) In-plane magnetic penetration depth 
in LSCO vs hole concentration S. 



D is related to the plasma frequency uj* by w 2 /8 = 
J °° du> Re cr(u>). In optical reflectivity measurements, 
the frequency-dependent conductivities can be derived 
from the reflectivity spectra. By integrating the spec- 
tral weight below 1.25eV, (w*) 2 was found to vanish 
linearly with the decrease of doping concentration in 
the low doping regime, and for optimally doped YBCO 
(fi-w*) 2 ~ A.heV along the a-axis (without the contribu- 
tion from the CuO chain) i.e. D ~ 145meVm& 

FigUni shows the RMFT and VMC results for Drude 
weight. Our results agree with those obtained with a fi- 
nite cut-off of the integration in Eq|^ to get rid of the 
contributions due to transitions from the ground state 
to the "upper Hubbard band"^ those results include 
the correction of order of O(Jft). The Drude weight 
increases almost linearly in the underdoped regime. 
Around the optimal doping, our results of the Drude 
weight is about QOmeV , in agreement with the optical 
reflectivity experimental data given in RefpS^. 



A. Drude Weight 

The Drude weight or the superfluid stiffness D is a 
measurement of superfluid condensation. In the linear- 
response theory, within the approximation that in the t-J 
model all optical spectral weights are condensed to zero 
energy, the Drude weight D can be given by^ 

(21) 



B. The antinodal quasiparticle energy E(n, 0) 

Around the antinodal point (7r,0), the quasiparticle 
dispersion becomes flat. This flat band has been stud- 
ied intensively by experimentsi 18 i 31 i 67 i 68 i 69 In the RMFT 
study, the energy of quasiparticle excitations at (jt,0) is 
given by: 

E(w, 0) = V(-4s**' -£) 2 +4A 2 . (22) 

Figllll shows the doping dependence of E(tv,0) ob- 
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FIG. 8: (Color online) In-plane magnetic penetration depth FIG. 9: (Color online) In-plane magnetic penetration depth 
in YBCO vs hole concentration 8. in Bi2Sr2Ca\- x Y x Cu20s+s vs hole concentration 8. 



tained in our calculation compared with the experi- 
mental resultsiSiSi. The experimental results obtained 
by AIPES (angle-integrated photoemission spectroscopy) 
and ARPES agree well with each other. In LSCO, the 
energy position of the flat band lies about 200 ~ 300meV 
below the Fermi energy for lightly underdoped state, and 
is shifted up to the Fermi level quickly with increasing 
doping, finally crosses the Fermi level at optimal doping 
or slightly overdoping. In Bi-2212, two branches of fiat 
bands (bonding band and anti-bonding band) were ob- 
served due to the bilayer splitting. They are determined 
by the low and high binding energies of the peak-dip- 
hump character. The average E(w, 0) of the two bands 67 
is shown in Fig lllT b). The bonding band has the same 
doping dependence as that in LSCO. The anti-bonding 
band lies much higher than the bonding band. 
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FIG. 10: (Color online) Drude weight D vs hole concentration 
5. The results from the variational theory given in Ref^ are 
for the Hubbard model with a finite cut-off to get rid of the 
contribution due to transitions from the ground state to the 
"upper Hubbard band" . In optical reflectivity measurements, 
the Drude weight is proportional to (<^p) 2 , for optimally doped 
YBCO (hujp) 2 ~ 4.56^ along the a-axis i.e. D ~ U5meV<^ 



In our calculation, without taking the bilayer coupling 
into consideration, our theoretical calculation of E(tt, 0) 
in Bi-2212 should correspond to the average E(ir, 0). 
Similar doping dependence is observed experimentally. 
E(ir, 0) approaches to the Fermi level with increasing 
doping, but does not get very close to the Fermi level 
even in the overdoped regime in contrary to LSCO. In 
VMC and PL1 simulation, the value of E(ir, 0) is much 
more closer to E(tv, 0) of Bi-2212 but is much larger than 
in LSCO. 



C. Chemical potential shift 

FigEl shows the electron chemical potential shift ju, 
compared with experimental data for LSCQii and Bi- 
22132a. In the RMFT, fj, is given by 



fx = H + 



i d (if Q 

N 85 



(23) 



The experimental data were deduced from the shifts of 
photoemission and inverse-photoemission spectra of the 
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FIG. 11: (Color online) The quasiparticle energy E(tv,0) vs 
hole concentration 5. 
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FIG. 13: (Color online) Nodal quasiparticle weight Z vs hole 
concentration S. Z = 1/(1+ A) was estimated in refi-, where A 
is the coupling constant estimated from the spectra function 
in ARPES 74 . The results from VMC simulation— is presented 
for a complementary comparison. 



particle coherent peak at the gap-nodes^ or from the for- 
mula Z = 1/(1 + A) where the coupling constant A can 
be extracted from the real part of self-energy Re E(fe, ui) 
of the spectral function^. In the RMFT analysis, the 
nodal quasiparticle spectral weight is equal to the renor- 
malized factor of the hopping term g t . Our theoretical 
results for the doping dependence of Z agree well with 
the experimental results. The nodal quasiparticle spec- 
tral weight grows almost linearly in the whole doping 
region as shown. The results from VMC™ are presented 
for a complementary comparison. 



FIG. 12: (Color online) Chemical potential shift Jl vs hole 
concentration 8. The experimental data were deduced from 
the shifts of photoemission and inverse-photoemission spectra 
of the core states of LSCO^i and Bi-2213ffl. 

core states of LSCO and Bi-2212. In LSCO the chem- 
ical potential shift Ji was found to be pinned close to 
zero energy in underdoped regime>Z2i£iiZ£ In Bi-2212, the 
chemical potential shift is not pinned at zero energy and 
shows a more rigid-band-like behavior. 

In our calculations, the chemical potential shift agrees 
qualitatively with the experimental data. It is also con- 
sistent with the result obtained by the exact diagonal- 
ization of the t-t'-J models. Furthermore, the shift is 
found to be larger in Bi-2212 than in LSCO in the entire 
hole doping range, in agreement with the experiments. 



D. Quasiparticle spectral weight 

FigEDshows the nodal quasiparticle spectral weight Z . 
In ARPES experiments the quasiparticle spectral weight 
Z can be deduced from the spectral weight of the quasi- 



E. Superconducting gap 

Experimentally the maximal superconducting gap A m 
can be measured by the thermal conductivity, ARPES, 
or other techniques. For example, from the thermal 
conductivity^, A m can be determined by assuming 
A m = Ukpv^ll with "universal" Fermi velocity vp2& and 
"weakly" doping dependent kp2&. In ARPES, the mid- 
point shift of the leading edge of the quasiparticle spec- 
tral at (%, 0) is approximately equal to A m . One can also 
determine A m by fitting the gap dispersion on the Fermi 
surface with the formula A(</3) = A m cos2(/3, where <f> is 
the Fermi surface angle Q&ZS/ZLZIL 

In Fig^l (a,c), our theoretical results of A m = 
A(cos k x — cos fcj/)|(7r,o) are shown and compared with the 
experimental data^*£i*£i*Z£. The doping dependence of 
A m agrees with the experiments, but the absolute val- 
ues are about twice larger than the experimental ones in 
YBCO (Bi-2212). 

FigEl (bid) compare the value of A$c = 3fA with 
the BCS gap Abcs — 2.14kBT c obtained by assuming 
T c = T c maj£ (l - 82.6(5 - 0.16) 2 ), (T c max = 95K for Bi- 
2212, 35K for LSCO)Z2. A SC and A BC S are roughly 
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FIG. 14: (Color online) Panels (a,c): The maximal supercon- 
ducting gap A m vs hole concentration 8. Panels (b,d): com- 
parison of the superconducting order parameter Asc = gtA 
in our RMFT calculation with the gap Abcs = 2.14fcsT c , 
where T c is estimated from T c = T c max (l - 82.6(5 - 0.16) 2 ) 79 
(T c max = 35K for LSCO and T™ ax = 95K for YBCO (Bi- 
2212)). 



trum of experimental data of low energy physical prop- 
erties in cuprates. In our theoretical calculations with 
both RMFT and VMC, the only parameters are the spin 
coupling and the hopping integrals of electrons on the 
CuOi plane, which are known quite accurately. We have 
found a qualitatively good agreement between the theory 
and the experiments on almost all the quantities we have 
studied, including the specific heat, the thermal conduc- 
tivity, the in-plane magnetic penetration depth, and the 
antinodal quasiparticle energy, the Drude weight, and 
the superconducting gap. The agreement on the doping 
dependences of these properties is remarkable except in 
the heavy overdoped regime. The major discrepancy is 
on the absolute values of some quantities, which may be 
attributed to the large value of Vi estimated in the the- 
ory. The comparison would be quite satisfactory quanti- 
tatively if one had used a theoretical value of V2 by 2 ~ 4 
times smaller, which indicates a possibility of overesti- 
mate of the gap. (It is known that the gap estimated 
in the VMC calculation is overestimated by a factor of 
2 or more). This discrepancy could also be due to the 
simplification of the model Hamiltonian or the approxi- 
mate wavefunction. More investigation will be needed to 
address these issues. 



proportional to each other. 
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